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Abstract 

We consider a multidimensional monostable reaction-diffusion equation whose nonlinear¬ 
ity involves periodic heterogeneity. This serves as a model of invasion for a population facing 
spatial heterogeneities. As a rescaling parameter tends to zero, we prove the convergence to 
a limit interface, whose motion is governed by the minimal speed (in each direction) of the 
underlying pulsating fronts. This dependance of the speed on the (moving) normal direction 
is in contrast with the homogeneous case and makes the analysis quite involved. Key ingredi¬ 
ents are the recent improvement [4] of the well-known spreading properties [32], [9], and the 
solution of a Hamilton-Jacobi equation. 
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1 Introduction 

We consider the Cauchy problem 


f dtu^ = eAu'^ + -/ f —: u^'] in (0, oo) x 

{P‘^) I £ \£ J 

[u^(0,x) = g(x) in 

where u will typically denotes a population density, and the nonlinearity f{x,u) is periodic in 
X € and of the monostable type. The parameter e > 0 measures the thickness of the diffuse 
interfacial layer, which will account for the invasion front of the population. Our goal is to study 
the asymptotic behavior — or the singular limit, or the sharp interface limit — of (P^) as e —>• 0. 

The reaction-diffusion equation in problem (P'^) arises from the hyperbolic space-time rescaling 
u^{t,x) := u (|, I) of the heterogeneous equation 

dtu = Au + f{x,u). (1) 

Let us emphasize that the understanding of the long time behavior of (1) is not equivalent to that 
of the sharp interface limit of (P^). Roughly speaking, the former one deals with the stabilization 
of the interface into a predetermined shape after a long time, whereas the latter one keeps the 
memory of the shape of the initial data. In other words, the singular limit procedure describes 
some transient states, during which geometry is quite relevant. 

Let us now state the assumptions on the nonlinearity f{x,u). Let Li,...,Lf^ be given positive 
constants. A function h : —>■ M is said to be periodic if 

h{xi,...,Xk + Pfc, ■■.,xn) = h{xi, ...,xn), 

for all 1 < fc < TV, all (xi, ...,Xn) S ■ In such case, (0,Pi) x • • • x is called the cell of 

periodicity. Through this work, we assume that 

for all u G R+, f{-,u) : R^ —>• K is periodic. (2) 

Our second main assumption on the nonlinearity / is the following. 

Assumption 1.1 (Monostable nonlinearity). The function f : R^ x ]R_|_ -^R is of class C^’°‘ in 
{x,u) and in u, and nonnegative on R^ x [0,1]. Concerning the steady states of the periodic 
equation (1), we assume that 

(i) the constants 0 and 1 are steady states (that is, /(•, 0) = /(•, 1) = 0 in R^); 

(ii) Vu G (0,1), 3a; G R^, f{x, u) > 0. 

{iii) there exists some p > 0 such that f{x,u) is nonincreasing with respect to u in the set 
R^ X (1-p, 1]. 

Notice that, if 0 < p{x) < 1 is a periodic stationary state, then p = 0 or p = 1. Indeed, 
since f{x,p) > 0, the strong maximum principle enforces p to be identically equal to its minimum, 
thus constant and, by (m), the constant has to be 0 or 1. Hence, under the above hypotheses, 
equation (1) is often referred to as the monostable equation. Typical examples are of the form 
f{x,u) = p{x)f{u), where p{x) is positive and periodic, and / is a homogeneous nonlinearity 
possibly of the following types: fi{u) = u{l — u) (Fisher-KPP), f 2 {u) = u^{l — u) with r > 1 (weak 
Allee effect), fsiu) = e“^/"(l — u) (Arrhenius nonlinearity), or f 4 (u) = M(e^““ — 1) (Nicholson’s 
blowflies equation). 

The monostable problem (1) arises in various fields of physics and the life sciences, and es¬ 
pecially in population dynamics models where propagation phenomena are involved. Indeed, a 
particular feature of this equation is the formation of traveling fronts, that is particular solutions 
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describing the transition at a constant speed from one stationary solution to another one. Such so¬ 
lutions have proved in numerous situations their utility in describing the dynamics of a population 
modelled by a reaction-diffusion equation. 

Equation (1) is a heterogeneous version of the reaction-diffusion equation 

dtU = Au + f{u), (3) 


with / of the monostable type. Among monostable nonlinearities, one can distinguish the ones 
satisfying the Fisher-KPP assumption, namely u i—^ is maximal at 0, the most famous example 

f{u) = fi{u) = u(l —m) being introduced by Fisher [14] and Kolmogorov, Petrovsky and Piskunov 
[24] to model the spreading of advantageous genetic features in a population. The KPP assumption 
means that the growth is only slowed down by the intra-specific competition, so that the growth per 
capita is maximal at small densities. Due for instance to the lack of genetic diversity at low density, 
this assumption may be unrealistic. To take into account such a weak Alice effect, one may use 
the growth function f{u) = / 2 (m) = u’'(l — u), r > 1. The nonlinearity f{u) = f 4 {u) = — 1) 

is commonly used [19] to explain oscillations of a population of Australian sheep blowflies, Lucilia 
Cuprina, described by Nicholson [29]. Let us notice that our work stands in the class of monostable 
nonlinearities, and therefore covers all these examples coming from population dynamics models, 
and the Arrhenius case f{u) = fsiu) = e“^/“(l — u) which comes from combustion models. 

Nevertheless, the environment is rarely homogeneous and may depend in a non trivial way 
on the position in space (patches, periodic media, or more general heterogeneity...), so that one 
should take into account heterogeneities. We refer to the seminal book of Shigesada and Kawasaki 
[30], and the enlightening introduction in [10] where the reader can find very precise and various 
references. For example such heterogeneities are very pertinent in some epidemiology models, 
where different treatments (antibiotics or insecticides) are tested, aiming at finding an optimal 
combination. 

In a periodic framework, traveling fronts in the homogeneous equation (3) are replaced by the 
so-called pulsating traveling fronts in the periodic equation (1) (see below for details). As far as 
the rescaled equation in (P*^) is concerned, fronts become sharper as e —>■ 0, and we therefore 
have to deal with the so-called interfaces. Also, as explained above, the singular limit analysis of 
(1) describes a transient state where the geometry of the initial habitat of the population is an 
insightful information. 

In this paper, we aim at looking at the way those interfaces are generated and propagate, 
hence providing some accurate connection between the behavior of solutions u^{t,x) in the fast 
reaction and low diffusion regime and some free boundary problem. One of the originality of 
this work is that we allow the equation to be spatially heterogeneous, which as recalled above 
is essential in realistic biological models. More precisely, we restrict ourselves to the spatially 
periodic case, which provides insightful information on the role and influence of the heterogeneity 
on the propagation, as well as a slightly more common mathematical framework. 

We will describe in subsection 2.1 what is known as far as front-like solutions of (1) are 
concerned. In particular, we will see that the outcome of the heterogeneity is some new dynamics, 
that do not appear in the homogeneous case, where the speed of the propagation depends on 
its direction. This feature is the origin of new technical difficulties when retrieving the interface 
motion. 


As far as initial data g{x) appearing in {P^) are concerned, we make the following hypotheses. 


Assumption 1.2 (Structure of initial data). fi) Let Dq he a nonempty, open and hounded set 
o/R^. Let g : Do —>■ [0,1) he a map of the class on Dp, positive on Dp and such that 
g{x) = 0 for all x G SDp. Define the map g : M. hy 


9{x) 


g{x) ifxGflo 
0 if x ^ flo . 


(ii) We assume that Dp is convex and has a smooth boundary Tp := dLlo- 
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Notice that the assumption g(x) < 1 becomes unnecessary if one assumes further that there 
is no steady state for (1) above 1. Also, rather than compactly supported initial data, one may 
allow g{x) to have tails that are “consistent” with those of the pulsating fronts (see [2] for the 
homogeneous case with “tails”). For the sake of simplicity, we do not consider here such cases. 
The convexity assumption (ii) will allow to describe explicitly the limit interface (obtained via a 
Hamilton-Jacobi approach) in Proposition 2.5 and then to use a family of planar supersolutions 
in Section 7. 

Before stating our results, let us now comment on related works. First, there is a large literature 
on the singular limit of (generalizations of) 

= eAu® + -f (x, . (4) 

e 

Observe that (4) arises after a hyperbolic rescaling of 

dtu = Am + /(ex, u), 

whereas Problem [P’^) under consideration follows from (1). First results are due to Freidlin [15, 16] 
using probabilistic methods. Later, Evans and Souganidis [13] used PDE technics, Hamilton-Jacobi 
framework to be more precise, to study (4). In this context, we also refer to [6], [7] and, for an 
overview, to [31]. Let us also mention the related work [27] which is linked with homogenization 
processes [25]. As far as (generalizations of) the considered problem (P^) is concerned, we refer to 
[26, Section 9] where Hamilton-Jacobi and homogenization technics are combined. Nevertheless, 
notice that all these results hold under the KPP assumption, that is f{x,u) < fu{x,0)u, whereas 
we stand in the larger class of monostable nonlinearities. 

In the homogeneous case f{x,u) = /(m), the sharp interface limit of (4) has been recently 
revisited using specific reaction-diffusion tools, such as the comparison principle and traveling 
wave solutions, which allows to capture accurate convergence rates [I, 2]. Hence, the introduction 
of a delay effect has been handled in [3], via such methods. 

Our analysis of the introduction of heterogeneity in {P‘^) stands mainly in this latter framework. 
It relies on accurate “local” subsolutions combined with improved spreading speeds properties [4], 
and on a family of planar supersolutions whose envelop solves the limit Hamilton-Jacobi equation. 


2 Some known results 

Before stating our main results in Section 3, we need to say a few words on monostable pulsating 
fronts and spreading speeds (in subsection 2.1), and on the limit free boundary problem (P^j) (in 
subsection 2.2), which is expected to describe the motion of the transition layers of the solutions 
u^{t, x) of (P^), as e: —>■ 0. 

2.1 Monostable pulsating fronts and spreading properties 

The definition of the so-called pulsating traveling wave was introduced by Xin [33] in the framework 
of flame propagation. It is the natural extension, in the periodic framework, of classical traveling 
waves. Due to the interest of taking into account the role of the heterogeneity of the medium on the 
propagation of solutions, a lot of attention was later drawn on this subject. As far as monostable 
pulsating fronts are concerned, we refer to the seminal works of Weinberger [32], Berestycki and 
Hamel [9]. Let us also mention [11], [20], [21], [28] for related results. 

For the sake of completeness, let us first recall the definition of a pulsating traveling wave for 
the monostable equation (1), as stated in [9]. 

Definition 2.1 (Pulsating traveling wave). A pulsating traveling wave solution, with speed c > 0 
in the direction n £ is an entire solution u{t,x) — t £ M, x £ — of (1) satisfying 

TT / k-n 

Vfc£j|LiZ, u(t,x)=uit-\ - ,x + k 
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for any t G R and x G , along with the asymptotics 

u{—oo, •) = 0 < u{-, •) < u{+oo, •) = 1, 

where the convergences in ±cxd are understood to hold locally uniformly in the space variable. 

One can easily check that, for any c > 0 and n G u{t,x) is a pulsating traveling wave 

with speed c in the direction n if and only if it can be written in the form 

u{t, x) = U{x ■ n — ct, x), 
where U{z,x) — z gR, x G — satisfies 

for all z G R, U{z, •) : R^ —>■ R is periodic, 

{7(— oo, •) = 1 <[/(•,•) < U{+oo, •) = 0 uniformly w.r.t. the space variable, 
along with the following equation 

{d^^ +A^)U+ 2V^d,U-n + cd^U + f{x,U) = 0 on R x R^. (5) 

We can now recall the result of [9], [32], on existence of pulsating traveling waves for the 
spatially periodic monostable equations: in any direction there is a minimal speed c*(n) > 0 
which allows existence. Precisely, the following holds. 

Theorem 2.2 (Monostable pulsating fronts, [9], [32]). Assume that f is of the spatially periodic 
monostable type, i.e. f satisfies (2) and Assumption 1.1. 

Then for any n G , there exists c*{n) > 0 such that traveling waves with speed c in the 
n-direction exist if and only if c> c*(n). Furthermore, any pulsating traveling wave is increasing 
in time. 


In the Fisher-KPP case the continuity of the velocity map n i—>■ c*{n), even if not explicitly 
stated, seems to follow from the characterization of c*{n) (see [32], [9]). In the more general 
monostable case, such a property was recently proved. 

Theorem 2.3 (Continuity of minimal speeds, [4]). The mapping n G >—>■ c*(n) is continuous. 


The introduction of these pulsating traveling waves was motivated by their expected role in 
describing the large time behavior of solutions of (1) for a large class of initial data. In this 
context, let us recall the seminal result of [32]: for any planar-like initial data in some direction n, 
the associated solution of (1) spreads in the n direction with speed c*(n). Actually, for our singular 
limit analysis, it turns out that we need the stronger property that this spreading is uniform with 
respect to the direction n. This was the purpose of our previous work [4]. 


Theorem 2.4 (Uniform spreading, [4]). Assume that f is of the spatially periodic monostable type, 
i.e. f satisfies (2) and Assumption 1.1. Let a family of nonnegative initial data ('Uo,n)nGS'^-i such 


that 


3C > 0, 
3/r > 0, 


Vn G 
3K > 0, 

inf inf 

nggN-i xGR'''' 


n>C uo^n{x) = 0, 
inf Mo,n(a;) > 

, x-n< — K 


We denote by (un)neS^-^ associated family of solutions of (1). 

Let a > 0 and p > 0 be given. Then, there exists r > 0 such that for all t > t, 


sup sup \1 - Un{t,x)\ < p, (6) 

x-n<{c* (n) — a)t 

sup sup Un{t,x) < p. (7) 

X'n'>{c* {n)-\-a)t 

Let us notice that, under suitable assumptions such as those in [9], [4], the above results 
are also available for more general spatially periodic and monostable equations which include 
heterogeneous diffusion and advection terms. We restrict ourselves to Problem (P^) to simplify 
the presentation, but our argument easily extends to such a framework. 
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2.2 On limit free boundary problems 


We recall that we aim at investigating the £ —>■ 0 limit of u^{t,x) the solution of (P^). Then the 
limit solution x) will be a step function, taking the value 1 on one side of a moving interface 
which we will denote by Ft, and 0 on the other side. This sharp interface, if smooth, obeys the 
law of motion 


(P°) 


Vn = c*{n) on Tt 


where Vn denotes the normal velocity of Tj in the exterior direction n, the unit outer normal of 
Ti at each point a: G Ft. Flere c*(n) denotes the minimal speed of the underlying monostable 
pulsating wave traveling in the n-direction. 

As we only know the mapping n i—>• c*{n) to be continuous, the smoothness of the interface, 
and hence the well posedness of (P°), is not guaranteed even for small positive times. 


A classical way to overcome the lack of smoothness is to define the limit interface via the level 
sets of the viscosity solution of a Flamilton-Jacobi problem 




+ |Vw|c* (^ 1 ^^ = 0 in(0,oo)x]R^ 
w(0,a:) = Wo{x) in 


Here wq : > R is any uniformly continuous function such that 

rio = {a: : woix) < 0}, Fq = {a; : wq{x) = 0}. (8) 

Thanks to the continuity of c*{n) with respect to n G namely Theorem 2.3, the Hamilton- 

Jacobi problem admits a unique viscosity solution w G C'((0,oo) x R-^), and 

fit := {a: : w{t,x) < 0}, Ft := {x : w{t,x) = 0} 

do not depend on the choice of wq as above (see Theorems 4.3.5 and 4.3.6 in [17]). As long as 
the solution of (P°) has a smooth solution, both motions coincide, which is why we still denote 
it by Ft. However, the Hamilton-Jacobi approach does not require smoothness as (P°) does, and 
therefore enables to define the zero level set Ft as the limit interface for all t > 0. 

The literature on this level set approach via viscosity solutions of Hamilton-Jacobi equations 
is rather large. The reader may consult [12] or the book of Giga [17] and the references therein. 

Thanks to the convexity of the initial set Hq, a so-called Hopf formula [23] is actually available 
and provides an explicit depiction of the motion, as stated in the following result. 

Proposition 2.5 (The limit interface explicitly). Let Assumption 1.2 {ii) hold. Let the limit 
interface Ft be defined via the Hamilton-Jacobi problem {P%j) as above. 

Then, for all time t>0, the set Ft is the zero level set of the convex function 


v{t, x) := max {x — y).ny — c*{ny)t, 
yeTo 

where Uy denotes the outward unit normal vector o/Fq at point y. In particular, for all time 
t > 0, the set Ft remains sharp, in the sense that it does not develop an interior, and the bounded 
domain fit delimited by Tt remains convex. 

Roughly speaking, this proposition means that the motion can be described by first looking 
at Fq as the envelop of some half-spaces, and by then letting each of those half-spaces move at 
the speed c*(n) corresponding to its normal direction. We refer to [5] where the Hopf formula 
was revisited in the context of viscosity solutions, and obtained using the more general theory of 
differential games. We propose a direct proof of Proposition 2.5 in subsection 4.1. 
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3 Main result 


We are now in the position to state our main result of convergence of to the interface motion 
defined via the level sets of solutions of {P^j)- Together with Proposition 2.5, the theorem below 
provides a precise depiction of the shape of solutions or, in other words, of the expansion of the 
habitat of the population. 

Theorem 3.1 (Convergence to a propagating interface). Let the nonlinearity f be of the spatially 
periodic monostable type, i.e. f satisfies (2) and Assumption 1.1, and let the initial data g in 
Problem (P^) satisfy Assumption 1.2. For any e > 0, let : [0, oo) x M.^ —^ R 6e the solution of 
(P^). Let Tt and Lit bs defined via the Hamilton-,Jacobi problem {P%j) as in subsection 2.2. 

Then, the following convergence results hold. 

(i) For any 0 < r < T < +oo and small /3 > 0, we have 

sup sup |1 — a;)| —>■ 0 as e —>• 0; 

T<t<T {x:d{t,x)< — ^} 

(a) For any 0 < T < +oo and small /3 > 0, we have 

sup sup a;)| —>■ 0 as e —>• 0. 

0<t<T {x-.d(t,x)>l3} 

Here d(t,-) denotes the signed distance to the set Pj, which is chosen to be negative in Lit o,nd 
positive in R^ \ (P* U fit). 

The rest of the paper is devoted to the proof of Theorem 3.1 and is organized as follows. 

We start, in Section 4, by some results on the motion of the limit interface which are crucial to 
our analysis of the parabolic problem (P"^), but are also of independent interest for the Hamilton- 
Jacobi problem (Pj^j). On the one hand, we prove Proposition 2.5, hence providing an explicit 
description of the limit interface. On the other hand, we approximate the motion defined via 
(P))j) by a smooth motion, which preserves all its essential geometric properties. 

To prove the control from below (i) of Theorem 3.1, we distinguish two regimes. First, we 
prove in Section 5 the emergence of transition layers for u^{t,x) in very small times. The prop¬ 
agation of the layers (from below) that occurs in later times is then studied in Section 6. The 
heterogeneity rises some technical difficulties since pulsating fronts depend non trivially on the 
direction of propagation. Roughly speaking, we construct “local” subsolutions and combine the 
uniform spreading properties of Theorem 2.4 with an iteration procedure. The construction of 
such subsolutions requires smoothness of the interface, which insures that the motion is locally 
governed by the planar dynamics of the rescaled equation (1). Hence, we actually apply the above 
procedure to the smooth approximated motion defined in Section 4. 

Last, in Section 7, to prove the control from above {if) of Theorem 3.1, we construct a family 
of planar supersolutions — whose envelop (4.1) coincides with the explicit characterization of 
Proposition 2.5 — and use again the uniform spreading properties of Theorem 2.4. 

4 Some results on the motion of the limit interface 

In this section, we are only concerned with the limit interface motion (P))j). We first prove 
the explicit description of Proposition 2.5, and then proceed to an approximation of the motion 
(Pf)j) by a smooth motion. As mentioned before, smoothness will play an essential role in the 
convergence of solutions of (P'^), and more specifically in Section 6. 

4.1 Characterization of the motion 

We begin by recalling that 

vit, x) := max {(x — y).ny — 
yeTo 


7 


c*{ny)t}, 


where Uy is the outward unit normal of Fq at point y. The zero level sets of x) are obtained by 
“intersecting all the half-planes arising from ?/ € Fq and propagating with speed c*{ny) in direction 
Uy”. We will prove that, at least for its level sets lying above some small —5 < 0, the function v is 
a viscosity solution of the Flamilton-Jacobi problem {P%j)- As the motion of interface is defined 
by the zero level set of the viscosity solution, this will be enough to infer that its zero level set 
defines the appropriate interface F 4 , that is Proposition 2.5. 

Remark 4.1. Write v(t,x) = maxygro where 

V' : (t, X, y) G (0, 00 ) X x Fq (x — y).ny — c*{ny)t, 

is continuous with respect to y G Fq, smooth and convex (since linear) with respect to t > 0 and 
X G For a given {t,x) G (0, 00 ) x let us denote by Y{t,x) the set o/y G Fq that maximize 
i({t, X, •), that is 

= {y G Fo : v(t,x) = 'ip{t,x,y)}. 

If, for a given (to,xo), the set Y(to,xo) reduces to a singleton yo then it follows from classical 
results of convex analysis (see [22, Corollary 4-4-5]) that v is differentiable at {to,xo), and 

dtv{to,xo) = dtip{to,xo,yo) = -c*{nyf), W„;v{tQ,xo) = W„ijj{to,xo,yo) = Uy^, 

so that V satisfies the Hamilton-Jacobi equation dtv |Vu|c* = 0 in the classical sense at 

(to,Xo). However, we have to deal with the case where Y{to,xo) is not a singleton. As we will see, 
this can be performed in the set {u > —J} for some small enough J > 0, and requires to cut-off 
the set {u < —J}. 

We prove the following proposition, of which Proposition 2.5 is an immediate corollary. 

Proposition 4.2 (A solution of the Hamilton-Jacobi problem). For any small enough J > 0, the 
function 

vs(t,x) := max(—(5; w(t, x)), 

is a (viscosity) solution of the equation of the limit problem {,P%j), that is 

-I-iVu^lc* = 0, m(0,oo)xR^, (9) 

and, by convexity, vs(0,x) is an admissible initial data for {P^j) the sense of ( 8 ). 

Proof. Recall that at time t = 0, Fq is a smooth hypersurface, and that the bounded set Hq 
delimited by Fg is convex. Hence, for J > 0 small enough, one can define a smooth hypersurface 
F--^ as 

Fq {x G : d{0,x) =-6} = {y - Suy : y G Fg}, 

where d(0, •) denotes the signed distance to Fg, which is negative in the bounded set Hg, and 
positive in R^ \ Hg. Notice also that, when x G Hg, we can write u(0, x) = — minygro dist{x, Hy), 
where Hy is the hyperplane going through y and with normal vector Uy. As a result, the convexity 
assumption yields 

Fq^ = {xG R^ : 'i;(0,x) =-J}. 

In particular, since the function v is convex with respect to x, the bounded set Hq delimited 
by Fj)"^ is still convex. Moreover, as we have chosen 6 small enough so that Fj)"^ is smooth, it is 
straightforward that the outward unit normal vector of Fq at y — Suy is also Uy. Therefore, by 
some slight abuse of notation, when y G F))"'^, Uy will denote the outward unit normal vector of 
Fq at point y. Then 

vs{t,x) = max{—(5, max{(x — y).ny — c*(ny)t}} 
yeFo 

= max{—(5, max{(x — (y — Sny)).ny — c*[ny)t — J}} 

yeTo 

= max{0, max {(x — y).ny — c*{ny)i\} — 5. 
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Therefore, vs{t,x) is a solution of {P^j) if and only if 

vs(t,x) := max{0, max {(a: — j/).ny — c*(ny)t}} 
yer-* 

is. For convenience, denote 

: (t, X, y) G (0, oo) x x F,^'^ >->• (x — y).ny — c*(ny)t, 

which is continuous with respect to y G smooth and linear with respect to t > 0 and x G R^, 
and introduce also 

ws{t^x):= max tjj^(t,x,y)^ 

so that vs(t,x) = iJiayi{0, ws{t,x)). 


Let us now prove that vs is a solution of (9). First, the null function and each function 
(t,x) I—>• il’^{t,x,y) solve (9) so that vs{t,x) — as a supremum of solutions — is a viscosity 
subsolution of (9). 

To prove that vsit, x) is also a supersolution, let be a smooth test function such that vs — ip 
has a zero local minimum at some point (to,xo) G (0,oo) x R'^. We need to prove that 


dtip{to,xo) + \\7(p{to,Xo)\c* 


/ Vtp(fo,xo) \ 
V|V(^(to,a:o)|/ 


> 0 . 


( 10 ) 


If ws{to,Xo) < 0, then = 0 in a neighborhood of {to,Xo) and (10) is clear. Let us now assume 
0 < ws{to,Xo) = vs{tQ,xo). Since vs — ip has a zero local minimum at (to,xo), the time-space 
gradient of (p at (to^xo) must belong to the time-space subdifferential of vs at {to,xo), which is 
given by 


dvs{to,xo) = 


dws{to,xo) 

Co {(0R,0Riv) U dws{to,xo)} 


if ws{to,xo) > 0 
if ws(to,xo) = 0, 


where Co A denotes the convex hull of the set A. It also follows from [22, Theorem 4.4.2] that 


dws{to,Xo) = Co {{dttlj\to,Xo,y),'^x'flj\to,xo,y)) = {-c*{ny),ny) £ R x R^ : y G F(to,a;o)}, 
where Y(to, xq) is the set of y £ Fq that maximize xq, •)• Hence, in any case, one can write 


p p 

dt(p(to,xo) = ^ -\iC*(ni), V(p(to,xo) = ^ \ni, 

i=l i=l 


for some yi,...,yp in Y(to,xo), and rii the outward unit normal of Fq at point yt, and some 
nonnegative Ai,...,Ap such that — f- Therefore our goal (10) is recast as 


* f A ELiA.c*(n0 

'vIELiA.n.i; - |ELiA,n.| 


Let us define 


no := 




£ § 


Af-l 


and pick a ?/o G Lg such that Uy^ = no- Note that such a yo necessarily exists from the smoothness 
of the bounded hypersurface F]]"^. One must then have 


tjj^(tQ,xo,yo) = (xq - yo).no - c*(no)to < ws(to,xo). 
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so that 


c*{no)to > {xo - yo).no - wsito^xo) 

LLi ^*(2^0 - yo)-?T- 


> 


lELi A.n,| 
LLi A*(a;o - yi)-'ni 
IELiA.n.| 


- Ws{to,Xo) 


- Ws{tQ,Xo). 


Here we used the convexity of , so that (yi — yo).ni > 0 for all 1 < i < p. Next, as each 
belongs to Y{to,xo), we have ws(to,xo) = (xq — yi).ni — c*{ni)to, so that 


c*{no)to > 


SLi AiC*K) 

IELiA.n.| 


to + Ws(to,Xo) 


ELiA. 

J2i=l Ai''T'i 


- 1 


ELi Azc*K) 

- IELiA.n.| 


since ws(to,xo) > 0 (notice that this is where it fails if no cut-off is performed). This proves (11) 
and concludes the proof of Proposition 4.2. □ 


4.2 Regularization of the motion 

We now construct, by the vanishing viscosity method, a smooth hypersurface r“ which approx¬ 
imates the interface Pj as a —>■ 0. Moreover, the motion of this smooth hypersurface is always 
“slower” than that of the original interface Pj: this will allow us, in Section 6, to construct 
subsolutions of (P*^) which fully cover the bounded set delimited by P“. 

Proposition 4.3 (Approximated smooth motion). Fix oq > 0 small enough and, for any 0 < 
a < ao, let P“ : M. be a smooth function such that 

0 < P“(p) < |p|(c*(p/|p|) - a), for all p G 

and, as a ^ 0, 

F°‘{p) —>• |p|c*(p/|p|), locally uniformly in . 

Let Vg (x) be a smooth and strictly convex function such that 

II Vuq Iloo < 1, -^5(0, •) + a < '*^0 < «<5(0, •) + 2a, 

where vs is the explicit viscosity solution of {Pfjj) with initial data vs{0,x), as defined in Propo¬ 
sition 4-2. 

Then, the solution v°‘ of the parabolic equation 

f dtv°‘ P“(Vu“) — aAv°‘ =0 in (0, oo) x 
|u“(0, x) = Uq (x) m R^, 

is smooth, convex w.r.t. space, and converges locally uniformly to vs as a ^ 0. In particular, for 
any T > 0 and up to reducing a, the zero level set P“ := {x G R'^ : v°‘{t,x) = 0} is a smooth 
hypersurface for any 0 < t < T, and is such that 

sup c??i(P“,Pt) —>■ 0 as a —>■ 0, (12) 

0<t<T 

where du{A,B) := maxjsupjjg^ dist(a, P), sup^ dist{b, A)} denotes the Hausdorff distance be¬ 
tween two compact sets A and B. Last, satisfies 

-f |Vx“| ^c* ^ 1^^^ - a^ > 0 m(0,oo)xR^. (13) 
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Proof. One can differentiate (in any direction) the parabolic equation satisfied by v°‘ and, using 
llVuJfIloo < 1 for any 0 < a < oq, deduce from the comparison principle that 

II Vu“(t, •)||oo < 1, for all 0 < a < ao,t > 0. 

In other words, the family ■))o<a<ao,t>o is uniformly Lipschitz-continuous. As confirmed by 

[18], the proof of Theorem 4.6.3 in [17] still applies thanks to the above estimate, even though the 
solutions we consider are unbounded. Therefore, one can conclude that the family of functions v°‘ 
converges locally uniformly to the unique viscosity solution of {P%j) with initial datum U 5 ( 0 ,x), 
namely vs- 

We now proceed by noting that, for each 0 < a < ooj the smoothness of follows from 
standard parabolic estimates. One can then differentiate the parabolic equation twice in any 
given direction e G and deduce from the comparison principle (recall that Vq is convex) that 

u“(t, •) is convex for any positive time. In particular, we have Av°‘(t,x) > 0 for all t > 0 and 
X € R^, which proves (13). 

Let us now turn to the convergence of the zero level set r“ of to Tj. The proof again follows 
the steps of [17] (see the proof of Theorem 4.6.4 in the particular case of geometric motions). We 
fix any /3 > 0 and T > 0 and show that, for small enough a, supg<(<jn d«(r?,rt) </3. By (13), 
we get that v°‘{t,x) > vs{t,x) for alH > 0 and x € in particular, r“ C ilt for all t > 0. Let 
now i? > 0 be large enough so that for all 0 < t < T the inclusion lit C Bn holds, where Bn 
denotes the ball of radius R and centered at the origin. By the locally uniform convergence, it is 
clear that for any small enough a and a; S fit such that d{t,x) < —/3 (recall that denotes 

the signed distance to Ft), then v°‘{t,x) < 0. The convergence (12) easily follows. 

Let us again fix T > 0 and now prove that, for small enough a, the zero level set r“ is a smooth 
hypersurface on the time interval [0,T]. Note that, for any Q <tQ <T and xq G r“^, one has that 
|Vu“(to,a:o)| ^ 0 provided a is small enough. Otherwise, it would follow from the convexity of 
•) that u“(to, •) > 0 in : a contradiction with the fact that it approaches vs {to, •) locally 
uniformly. Then, as |Vu“(to,a^o)| 7^ Oj one can apply the implicit function theorem and obtain 
the smoothness of . □ 

5 Rapid emergence of the layers from below 

In this section we prove that, as e —>■ 0, the solution u‘^{t, x) of {P'^) is very close to 1 in fig after a 
very short time. The proof relies on the spreading properties of solutions of (1) with large enough 
compact support at initial time [32]. Precisely, the following holds. 

Proposition 5.1 (Emergence of the layers from below). Let the nonlinearity f be of the spatially 
periodic monostable type, i.e. f satisfies (2) and Assumption 1.1. Let the initial data g in Problem 
(P^) satisfy Assumption 1.2. 

Then, for any small rj > 0 and small a > 0, there is a time ta > 0 such that the following 
holds: there is eg > 0 such that, for all e G (0,£o), 

X G flo) dist(x,dLlo) > a l — g< u^{tae,x) < 1. (14) 

Proof. Since 1 solves the reaction-diffusion equation in (P'^) and since ^^(0, •) = g{-) < 1, the 
comparison principle implies u‘^{t,x) < 1, which proves the upper bound in (14). We next prove 
the lower bound. 

We begin by recalling the following result on the spreading of solutions with initial compact 
support [32, Theorem 2.3]: for any a G (0,1), there is R^ > 0 large enough so that the solution 
V of (1) with initial datum vq = ctxbr,, converges locally uniformly to 1 as t —>■ -l-oo. Here, x 
denotes the characteristic function and Bn the ball of radius R and centered at the origin. Note 
that Weinberger’s result [32] also provides a positive spreading speed in any direction; however, it 
is not required to prove Proposition 5.1. 
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Let us now fix some 77 > 0 and a > 0. From Assumption 1.2 on the initial data g, there is 
CTi G (0,1) such that, for all £ > 0, 

X G flo, dist{x,dV,o) > a => u^{0,x) = g{x) > ai. (15) 

We can now let to, > 0 be such that the solution u of (1) with initial datum vg = satishes 

v{ta,x)>l-g, VxGSafl^^. (16) 

We assume without loss of generality that > 2'/N max^ Li. 

Let us now fix x* G Oq such that dist{x*, dflg) > a. We are going to prove 

u^{tae,x*) > 1 — g, (17) 

for £ G (0,£o), where £0 > 0 has to be independent on the point x* chosen as above. We let 

xg G dflg such that dist{x*, dflg) = |a;* — a;o|. Since Rai > 2y/N max^Li, there exists k* = 

(fc* g,.., fc)(f g) G such that 

X* - 2R^,e ,^°~^ , G ek*L + Br^,s, (18) 

\Xg-X*\ 

where we denote k^L := (fc* ^L^r). Also, provided a and £0 > 0 are small enough 

depending only on 0 < maxygro l{y) < +00 with 7(2/) the mean curvature (positive by convexity) 
of Fg at point y, we have, for all £ G (0, £0), 

£k*L € Q,g and dist(ek*L,dLtg) > a + Ra^^e. (19) 

Observe that if x ^ ek*L + then uq ^^at if 

X G £k*L + Bsr^_^ then (19) implies that x G Op and dist{x,d^}g) > a. Hence, it follows from (15) 
that 

g{x) > V ^0, -—for all x G R^. 

Since v{^, ^ ek^L -^ solves the parabolic equation in (P'^), the comparison principle implies in 
particular that 

_, / X* — ek*L 

U [ta£, X ) > V [ta, -^— 

\ e 

In view of (16) and (18), the above estimate implies (17). The proposition is proved. □ 

The above argument also shows that, roughly speaking, the solution of (P^) may only expand, 
which is rather natural from the dynamics of the monostable equation. Precisely the following 
holds. 

Lemma 5.2 (Expansion). Let g > 0 be given. Let {Llt)g<t<T be a family of bounded and convex 
domains with smooth boundaries F; := dLlt. Then, for any a G (0,1) there is a time to- > 0 such 
that the following holds: there is eg > 0 — depending only on 0 < maxo<t<T max^^p^ lt{y) < +00 
with jtiy) the mean curvature ofTt at point y — sueh that, for any 0 <tg <T, any e G (0,£o), 
U^(tg,x)>(J, Vx G Hto => M^(t, x) > 1 — 77, Vx G iito, Vt > tp + 

6 Propagation of the layers from below 

We now begin the analysis of the motion of interface. In this section, we prove the lower estimate 
on the motion of level sets of the solutions n®(t,x), namely statement (i) of Theorem 3.1. 

To that purpose, we fix some times 0 < r < P, and a small /? > 0. We then let a > 0 be small 
enough so that the hypersurfaces (F“)p<t<7'+i, as defined in subsection 4.2, are smooth and such 
that 

sup d„(r“,ro<§. (20) 

0<t<T+l ^ 

We also denote, in this section, by the region enclosed by r“. 
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6.1 Lower estimates in small canisters 


We start by looking, for any fixed time to, at the “local motion” of the interface. By “local 
motion”, we mean that we will investigate the motion of the solution on small neighborhoods of 
any point of r“^. Precisely, the following holds. 

Lemma 6.1 (Lower estimates in small canisters). Let rj > 0 be given. Fix some time to £ (0,T), 
and assume that 

X £ n)), => x) > I — rj. (21) 

Then there are two positive constants Ai and A 2 , independent on to and e > 0 (provided it is 
small enough), sueh that 

u^{to + Aiy/e, x) > 1 — ry, 

for all X € D := Ua;j,grf C{xo), where C{xo) is the finite cylinder, or canister, made of the points x 
such that 

\{x — xo) ■ n\ < Ai (^c*(n) — y/e and |(a; — xq) • n'*‘| < ^ 2 ^. ( 22 ) 

Here n denotes the unit outer normal ofVf^ at point xo, and {x — xo) ■ n'^ denotes the orthogonal 
projection of x — xq on the hyperplane (Rn)-*-. 

Proof. First, let 7 > 0 be large enough so that, for all t £ [0, T], all y £ r“ with Uy the associated 
unit outer normal, we have the inclusion 


B^ (^y - C ilf, (23) 

where Br{z) denotes the open ball of center z, radius r. By convexity, it suffices to take 7 as the 
maximal curvature (in absolute value) of r“ in the time interval [ 0 ,T]. 

Let 77 > 0 and 0 < to < T’ be given. Let xo £ r“^ be given and n the associated unit outer 
normal. For the lemma to be proved notice that constants Ai and A 2 , that we need to determine, 
have to be independent on to, small e > 0 but also on xq and n. By assumption (21) and inclusion 
(23), we have 

Vx £ Bi (Xo - n ) , u^{to,x) > 1 — 77. 

^ V 7 / 

We fix a constant C > 2y/N max^ Li and, proceeding similarly as in Section 5, we can find 
some point ekeL := e(fci_eLi,..., k^^^LN), where kt^g £ Z for all 1 < 7 < IV, and such that 


Xo -£ skgL + Bce- 

7 


Then 


Vx £ Bi_q^ {skgL), u^{to,x) >1 — 77. 
This leads us to study the solution u(t,x) of (1) with initial datum 


(24) 

(25) 


Mo(x) := (1 - 77 ) X XB 1 (x), (26) 

where B^. denotes the open ball centered at the origin and of radius r. Note that this initial datum 
has compact support, so that Theorem 2.4 does not apply. In fact, the solution u(t,x) does not 
spread with speed c*(n) in the n-direction as t —> + 00 , but rather with some minimum of the 
over all n' £ However, as the radius of the initial support is very large, we can exhibit 

some transient dynamics where the solution does spread, in any direction n, with speed c*{n) the 
minimal speed of pulsating traveling waves. Let us make this sketch precise. 

We first note that, provided that e is small depending only on C and 7 , the finite cylinder 

Do := I X £ : lx • 77 ! < -2C and lx • n-*-! < 

7e 
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is a subset of thanks to Pythagoras’ theorem. In order to apply Theorem 2.4, which is 

•ye 

concerned with planar-shaped initial data, it is more convenient to consider a box-shaped initial 
support. With this in mind, we introduce (ni, ...,n 7 v-i) an orthonormalized basis of (Rn)-*-, and 
define the finite box 


Di :=< X G 


pAT 


li • n| <- 2C and VI < i < iV — 1, la; • nil < 

'ye 


C 


{N — l)7e 


which is a subset of Dq. 

We can now begin our investigation of the spreading of u, the solution of (1) with initial datum 
(26). By the parabolic comparison principle, we have 


u > u 


where u is the solution of ( 1 ) with initial datum 

Uo{x) := (1 - ??) X xdA^)- 

We let u{t,x-,n) denote the solution of (1) with initial datum 


uo{x-,n) := (1 - r?) X X{x-n<^-2C}(x)^ 

which is planar-shaped so that u{t,x;n) spreads in the direction n with speed c*(n). Precisely, 
recalling C > 2'/N max^Lj, we can find some point k^L := (fei^gLi,..., where ^ 

for all 1 < * < iV, and such that G Bc{keL). Then observe that 

Vo{x-n) := uo{x + hL-,n) > {1 - rj) x X{x-n<-3C}{x). 

We can now apply Theorem 2.4 with the family of functions in the right-hand side member above 
(which do not depend on e) as the family of initial data. Then, applying the comparison principle, 
we get that there exists r > 0 (which does not depend on e) such that 

inf inf v(t,x;n) > 1 — —, 

x-n<{c* {n) —^a)t 2 

where v{t,x;n) denotes the solution of (1) with initial datum vo(x]n). Then, since v{t,x\n) = 
u(t,x + keL]n) thanks to the spatial periodicity, the above estimate implies 

inf inf u{t,x;n) > 1 — (27) 

x-n<^—3C+{c''{n)—joi)t 2 


We emphasize that r > 0 can also be chosen independently of n € this is the exact purpose 

of our improvement of Weinberger’s spreading result [32], namely Theorem 2.4. 

We now estimate the difference w := u — u> 0, which satisfies dtw — Aw — g(t, x)w = 0, where 


g{t,x) 


f{x,u) - f{x,u) 
u — u 
duf{x,u) 


if w{t, x) ^ 0 , 
if w{t, x) = 0 . 


From Assumption 1.1, g(t,x) is uniformly bounded by some K which only depends on /. Then w 
satisfies 

dtw — Aw — Kw < 0. (28) 

As this parabolic equation is linear, we infer that w{t,x) < Wiit^x), where wq is the 

solution of (28) with initial datum 


wo(0,x) 


1-g if X-n <+2C, 
0 otherwise. 
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and the Wi’s, 1 < i < N — 1, are the solutions of (28) with initial data 


W2i-i(0,a;) = 


1 — 77 it X ■ n < ^ — 2C and x ■ rii > 
otherwise, 


C 


W2t{0,x) = 


1 — n it X ■ n < — - 2C and x ■ rij < — 

I — ■ye * — 

otherwise. 


{N — l)ye ’ 


C 

{N — l)ye ’ 


Note that, for any e € ^ and any positive constant M, (t,x) >—>■ Me ^VWt) 

supersolution of the linear equation (28). It therefore follows that 

Wo {t,x) < 

and, for any integer 1 < i < fV — 1 , 

-1 {t,x) < 

{t,x) < + 


W2i-1 


W2i I 


Then, we conclude that 


0 < (u — u) 


— u) f 


\\/s 


A 

v^’ 


2N-2 


E 


i=0 


E 

ve 


X < 


2 ’ 


(29) 


where 


^1 := T 


C 


4 y i^(Af-1)7’ 
for all X satisfying the two following inequalities: 


X ‘ Ti ^ -T 2C -f 2jt\ 

ye 


\x ■nA< 


C 


{N — 1)76 


- 2yli 


K 1 , / ?7 \ 1 ^ / 1 

e v^'^VdivJ" ye+^VV^ 

^ ^In 


(30) 


(^ "l 


\mJ 

2\ 


C 


0 ( 1 ), 


for 1 < * < — 1. The second inequality is in particular satisfied, for e > 0 small enough, if 

1 


\x ■ n^\ < 


C 


3\{N- l) 7 e ‘ i/e' 


(31) 


Combining the spreading property (27) of u and inequality (29), we conclude that 

u{te,x)> l-fj, te := E, (32) 

Ve 

for any x satisfying both inequalities (30) and (31), as well as 

X ■ n < - 3C + (c*(n) — 7 Q:)te. (33) 

ye 4 

We can now go back to our original problem (P^). Notice that both u(|, |) and u^{to+t, ekeL+ 
x) solve the equation in (P'^). Using Di C Pj__c and (25), we see that m( 0, f) < u^(to,ek^L + x) 

-ye ^ 

so that 

u^{to +1, ekeL + x) > u 
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where ekgL satisfies (24). Thus, we get 

u‘^{tQ + ete,x)>u(^e,- —(34) 

provided that satisfies (30), (31), (33) (so that (32) holds). Now, assume that x satisfies 

(22). Combining the first part of (22) and (24), we see that satisfies both (30) and (33). 

Combining the second part of (22), n • n-*- = 0 and (24), we see that satisfies (31). Hence, 

(34) holds true and is the desired conclusion that u'^{to + Axsje, x) >1 — rj. 

Note that, as announced, the constants Ai and A 2 defined above depend neither on to G (0, T), 
xq G r“^ and the associated unit outer normal n, nor on e > 0. The lemma is proved. □ 

Remark 6.2. Let us notice that Lemma 6.1 shares some ideas with the so-called consistency 
assumption (H4) of Barles and Souganidis [8]. Roughly speaking, their method consists in reducing 
the study of the sharp interface limit to compact and smooth shapes as well as to small times, that 
is to consistency. In a heterogeneous and bistable context, they then proved consistency under the 
additional assumption that the traveling wave (which in such case is unique) depends regularly on 
its direction. However, such a property is far from trivial, especially in the monostable case. We 
therefore adopt a different approach, relying on the uniform spreading properties proved in our 
earlier work [^], namely Theorem 2.4. 

6.2 Lower estimates for propagation of the layers 

We now complete our argument by combining an iteration method and Lemma 5.2. 

Proof of statement {i) of Theorem 3.1. We need to show that, for e > 0 small enough, we have 
u^^it, x) > 1 — r], for all r < t < T, for all x such that d{t, x) < —(3 (recall that d{t, •) denotes the 
signed distance function to Tj, negative in Hj). 

Recalling that Tq C LIq and a < d^(rQ,ro) < 2a (see Proposition 4.3), it follows from 
Proposition 5.1 that, for e > 0 small enough, assumption (21) of Lemma 6.1 is satisfied at time 
to = taS < T. As a result 

u^{ti,x) > 1 - r], ti := to-G Aiy/e, (35) 

for any x G D defined as in Lemma 6.1. Moreover, (35) also holds true if x G 0“^ in virtue of 
Lemma 5.2 (notice that the needed time to reach 1 — 77 in Lemma 5.2 is of order e), with r“, 
playing the roles of Pt, Lit. Therefore, it follows from the claim 

(36) 


(whose proof is postponed), that 


Vx G 0 “^ , u^{ti,x) >1 — 77. 

Proceeding by induction, we conclude that for all times 


tk ■= to + kAiy/e, 


up to some k such that T < < T + 1, we have 

u^{tk,x) > 1 — 77 for all X G 

In particular it follows from (20) that u^{tk,x) > 1 — 77 for any x such that d{tk,x) < —f3. 

It now only remains to consider intermediate times. Notice that, even though we stated in 
Lemma 5.2 that the solution may only expand, this is in fact only true when looking at 
interval of times of order larger than e. Therefore, the above inequality does not guarantee that 
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d{t, x) < —/? => x) > 1 — 77 in intervals of time [tk, tk + 0{£)\. To avoid this difficulty, we can 
nevertheless note that for all t € \tk,tk+i) with fc > 1, 

u^(t,x)>l —77 for all (37) 

Note that up to reducing e, we can assume that ti < t. Let now any t G and k > \ 

such that t G [tk,tk+i). Let also a; € fit be such that d{t,x) < —/3. Notice that it follows from 
Proposition 2.5 that there is C > 0 such that c?^(rt, rt^_J < C(t — tk-i) < 2AiC^. Recall also 
that a > 0 was chosen such that (20) holds, so that (iH(rt,r“) < |, for all r < t < r + 1. As a 
result 

c^w(rt,rtj_J < 2AiC-v/e + ^ < /3, 

for £ > 0 small enough. Since d{t,x) < —/3, this enforces x G ^y (37), we get that 

u^{t,x) >1 — 77 . This concludes the proof of the lower estimates on the motion of the layers of 
u^{t,x). □ 

Proof of claim (36). Recall that (see Proposition 4.3) r“ is the zero level set of •), where 

V7;“ \ \ 

To prove the claim (36), consider any x G \ and let us prove that x G D. First, there 
exists some Xq G P)), such that |a; — a;o| = dist(x,r“^) and, by convexity, such an xq is unique. 
Moreover, 

X - xo Vv°‘(to,xo) 

71 ^ - — - 

|a;-xo| |V7;“(to,a;o)| 

is, by construction, the unit outer normal of Tf^ at point xq (the first equality follows from the 
choice of a:o, and the second from the definition of Tf as the zero level set of 7;“(f, •)). 

In order to prove that x G D, only remains to check the inequality 

|(a; - Xo) • 77-1 = |x - xo| < Ai (^c*(n) - y/l. 

Note that, by convexity of 7 ;“, 

> v°(to,xo) + V?;“(to,a7o) • (x - xq). 


,9t7;“ + |V7;“| c 


and also that, thanks to the smoothness of 7 ;“, 

7;“(ti,x) - v°‘{to,x) > dtV°‘{to,x){ti - to) - K\ti - toP, 

where K is a positive constant (recall that a > 0 has been fixed). Since x G we have 
7;“(ti,x) < 0 and since xq G we have 7;“(to,a;o) =0. As a result, up to increasing K if 
necessary, 

0 > 7;“(ti,x) - 7;“(to,xo) 

> V7;“(to,a;o) • (x - xq) + dtV^{to,x){ti - to) - K\ti - toP 

> V7;“(to,a;o) • (x - xq) + {dtV°‘{to,xo) - K\x - xo|)(ti - to) - K\h - 

Using (38), we deduce that 

0 > |V7;“(to,a:o)| x |x-xo| - |V7;“(<o, a:o)| {c*{n) - a) {h - to) - K[\ti - to\ x |x-xo| + |ti -ffiP]- 
Recalling that V7;“ does not cancel on Tf, we can infer by compactness that 

p := inf inf |V7;“(t,x)| > 0. 

0<7<Ta:Grf 
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Recalling also that ti — to = it follows from the above that 


Ix-xol < 
< 


|V?;“(to,a;o)| 


|Vii“(to,a;o)| - KAi^/s 
{c*{n) - I) AiVe, 


(c* (n) — a)Ai yfe + 


KA\e 


p - KAx^fe 


provided e > 0 is small enough. As announced, x € D and the claim is proved. 


□ 


7 Control of the layers from above 

In this section, we prove the upper estimate on the motion of level sets of the solutions M^(t,x), 
namely statement (ii) of Theorem 3.1. 

To do so, we are going to construct a family of planar supersolutions (indexed by y € Tq) for 
(P®), whose envelop is close to the zero level sets of v{t, •), that is Tt in virtue of Proposition 2.5. 
Then, for the sake of clarity, rather than using the uniform upper spreading speed (7), we instead 
use some kind of uniform asymptotics of the monostable minimal waves — which is proved in [4] 
and actually implies (7). 

Lemma 7.1 (Uniform asymptotics for critical waves, [4]). Let u*{t, x;n) = U*{x-n — c*{n)t, x; n) 
be a family of increasing in time pulsating traveling waves of (1), with minimal speed c*(n) in 
each direction n £ , shifted so that U*(Q, 0;n) = ^. 

Then, the asymptotics U*{—oo,x;n) = 1, U*{oo,x;n) = 0 (which are uniform with respect to 
X £ ) are uniform with respect to n G 

Proof of statement {ii) of Theorem 3.1. Let 0 < T and a small /3 > 0 be given. For any n £ 
denote by U*{z,x;n) a monostable pulsating front with minimal speed c*(n) in the direction n, 
shifted so that t/*(0, 0;n) = 

Thanks to ||y||oo < 1 (see Assumption 1.2) and the above lemma, we can select some K > 0 
large enough so that 

U*{z, x; n) > llylloo, Vz < -AT, Vcc £ Vn £ (39) 


Then, for any y £ Tq and denoting again by Uy the outward unit normal of Tq at point y, we 
define 


u{t, x) := U* 



y).ny - c*{ny)t 

e 



From equation (5) for the traveling front, we deduce that u{t,x) solves the parabolic equation in 
(P®). We also have u^{0,x) = g{x) < u{0,x): indeed, for x ^ Uq we have g{x) = 0, whereas for 
a; £ Uq we have {x — y).ny < 0 by convexity and (39) gives the desired ordering. The comparison 
principle then implies u^{t,x) < u{t,x). As a result 


0 < u^{t,x) < inf U* 
veTo 



y).ny - c*{ny)t 
s 



(40) 


We recall that d{t, •) denotes the signed distance to the set Ft, which is chosen to be negative 
in Lit and positive in \ (Fj U 12*). Let us now prove that there is some 6 > 0 such that, for any 
t £ [0, T] and any x such that d{t, x) > (3, then 


3y £ Fo, {x - y).ny - c{ny)t>e(3. (41) 

Assume by contradiction that there are some sequences {tk)k>i, {xk)k>i as above such that 

Vy£ro, {xk - y).ny - c*{ny)tk < ^. 
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This enforces the sequence (a;fc)fc>i to be bounded so that, after extraction of some subsequences, 
we are equipped with some too G some Xod with d{tao,Xoo) > /? > 0, such that 

Vy G hoj (^oo y)'^y ^ (^y)too ^ 0- 

Thus v{too,Xoo) < 0, which contradicts d(too,a;oo) > /?■ 

Let us now choose any t G [0,T], any x such that d{t,x) > /3. In view of (41), we can select 
some yo G Tq such that (x — yo)-'ny„ — c*(ny„)t > Of}. Then, using (40) and the monotonicity of 
the pulsating traveling wave U*{z,x]n) with respect to its first variable, we get 

0 < u^{t,x) < U* (— — K,—;nyg] < sup sup U*(——K,X;n 

Thanks to Lemma 7.1, this implies supg<j<y ^'^P{x:d{t,x)>i 3 } t 0 as e —>• 0, which 

concludes the proof of Theorem 3.1. □ 
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